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We report a doping, magnetic field and low-temperature dependent study of the specific heat of 
the iron-arsenide Ba(Fei_a;Coa;)2As2 at under (x=0.045), optimal (x=0.08) and overdoped (x=0.103 
and 0.105) regimes. By subtracting the lattice specific heat the temperature and magnetic field de- 
pendence of the electronic specific heat has been studied. The temperature and field dependencies 
of the superconducting part of Cp exhibit similar behavior for all doping concentrations. The tem- 
perature variation of the electronic specific heat as well as its field dependence cannot be described 
by a single isotropic s-wave gap, pointing to a complex gap structure in the system. The lack of 
doping dependence indicates that the gap structure does not change significantly as a function of 
doping. We also observe a significant residual linear term of unknown origin in the specific heat 
of Ba(Fei_a;Co2,)2As2 which suggests that inhomogeneity may be an important factor in Co-doped 
BaFe2As2. 

PACS numbers: 74.20.Rp, 74.70.Dd, 74.62.Dh, 65.40.Ba 



The discovery of superconductivity in FeAs-based 
i?FeAsO[l|, (i?-rare earth) has opened a new era in su- 
perconductivity studies. Shortly after, other types of su- 
perconducting materials containing FeAs layers were dis- 
covered including: binary chalocogenides Fei+a;SeJp, 
so called "111" compounds LiFeAs or NaFeAs [1, [a|and 
122-systems ^Fe2As2 where A is an alkali element • 
Despite a large theoretical [l3| and experimental [ll| effort 
in the newly discovered Fe-based superconductors [l|, [2] , 
the nature of the superconductivity in these materials in- 
cluding the pairing mechanism and the symmetry of the 
order parameter remain unknown. Moreover, the experi- 
mental results reported so far are often contradictory, not 
only between various techniques, but also between differ- 
ent families. The large sample and doping dependence 
may favor scenarios where the low energy excitations, 
possibly nodal, strongly depend on the particular sam- 
ple being studied and the probe used to investigate them 
(e.g. RefUMl. 



Recently, much attention has been focused on the Co- 
doped BaFe2As2 familyQ due to the large single crys- 
tals which can be produced. They also appear to be 
more homogeneous than alternative dopings such as K- 
doped BaFe2As2[ll|. At optimal Co-doping (x = 0.08) 
an isotropic gap has been postulated by angle resolved 
photoemission spectroscopy (ARPES)[lal and scanning 
tunnel microscope (STM) JL7J measurements , while other 
experiments such as penetration depth [l8|. /uSR18|, 
NMRf2^, thermal conductivity [HHll], specific heat[2J, 



and Raman scattering.25,j point to an anisotropic gap 
scenario. Several of these measurements are consis- 
tent with the so-called s± model with a sign rever- 
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sal of the order parameter between different sheets of 
the Fermi surface [26l - [28| . Recently, it has been sug- 
gested by low-temperature thermal conductivity stud- 
ies, that in the Ba(Fei_a;Coa;)2As2 system the supercon- 
ducting gap evolves from uniformly large everywhere on 
the Fermi surface for low doping to having a very small 
value somewhere on the Fermi surface for high cobalt 
concentration[2l|. Raman scattering measurements have 
been made which support this conclusion [25I]. A simi- 
lar situation has been argued to exist in P doped FeAs 
compounds [2^. 

In this paper, we present results of our detailed studies 
of the specific heat of the Ba(Fei_a;Coa:)2As2 (x = 0.045, 
0.08, 0.103 and 0.105). By subtracting the lattice contri- 
bution, we extract the full electronic T-dependence for all 
compositions studied. A relatively large residual specific 
heat is related to the presence of the non-superconducting 
fraction in the samples. The temperature and field vari- 
ations of the superconducting part of Cp exhibit com- 
mon behavior for all doping concentrations indicating an 
anisotropic gap structure whose gross features are doping 
independent over the range investigated here. 

Single crystals were grown out of FeAs flux with the 
typical size of about 2x1.5x0.2 mm'^Q. The samples 
crystalize as well- formed plates with the [001] direction 
perpendicular to the plane of the crystals. The doping 
level was determined by microprobe analysis. The heat 
capacity was measured down to 400 mK and in magnetic 
fields up to 9 T using a thermal relaxation method im- 
plemented in a Quantum Design PPMS-9 device. All 
specific heat data measured in field were field cooled. 
Magnetic susceptibility have been taken in field cooled 
conditions with a field of 20 Oe applied parallel to the ab 
plane of the single crystals. Data on all samples of similar 
size and shape were normalized by a constant diamagne- 
tization factor which gave l/47r for x=0.08. 

In general, in the FeAs-based superconductors it is 
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challenging, due to the high iJc2, to obtain the nor- 
mal state electronic heat capacity in the superconduct- 
ing regime. In order to evaluate the electronic con- 
tribution of the specific heat, we have used a similar 
approach used previously for the optimal doped com- 
pound Ba(Feo.92Coo.o8)2As2 (see Reflsol. Isil ). We as- 
sume that the phonon part of the specific heat is inde- 
pendent of doping and we use the phonon specific heat 
obtained from the parent compound. BaFe2As2 shows 
a SDW transition at about 140 K. Recent inelastic neu- 
tron scattering experiments show that, in the ordered 
state, spin-wave excitations have a large gap of about 
10 meV (A « 116 K}M. Therefore, below 40 K, C,nag 
is almost negligible [3^ and we separate the lattice con- 
tribution to the specific heat of the parent compound 

as Cph = C^'^P'^^As^ - ^BaFe^As^rj. ^^^^^ ^BaFe^As^ jg 

the intercept of C/T of BaFe2As2. Thus, the 

electronic specific heat, at finite doping, is determined 

by (7g;(y)-S''(-P''=i-^C'°^)2^S2 _ Q^^rp-^Ba{Fei-^Co^)2As2 _ 

C{T)ph- Moreover, a small Schottky-like contribution of 
about 0.5 mJ/mol of the total specific heat at 0.5 K 
has been also subtracted from the data. The obtained 
temperature dependence of the specific heat for all sam- 
ples is shown in Fig.l. For the normal state specific heat 
below Tc we assume the form (C — Cph)lT = 7„ + hT. 
The hT term represents a small correction to the nor- 
mal state specific heat below Tc required for the samples 
not at optimal doping in order to conserve entropy be- 
tween the normal and superconducting states at T^. The 
so-derived normal state specific heat is shown by dashed 
lines in Fig.l. We cannot determine whether the nec- 
essary correction to the normal state specific heat im- 
plies either that our assumption of a doping independent 
phonon contribution is incorrect, that magnetic contribu- 
tions are non-negligible at some dopings, that a pseudo- 
gap is present, that quantum critical fluctuations exist, or 
that some combination of these effects are at play. How- 
ever, we emphasize that our conclusions are independent 
of the particular form of the normal state specific heat 
constructed to conserve entropy. 

As can be seen from the Fig.l, at very low temper- 
atures, for all compositions a significant residual spe- 
cific heat coefficient 70 is observed. It ranges from 
70 = 3.7 mJ/mol for x = 0.08 to 14.6 mJ/mol 
for X — 0.105. The values of the residual specific heat 
7o, normal state specific heat 7„ and the difference (7^- 
7o) as a function of Co concentration are presented in 
the inset of Fig.l. The smallest value of 70 is ob- 
served for the concentration close to the optimal Co dop- 
ing. The residual specific heat strongly increases when 
moving towards under or overdoped directions. Sim- 
ilar behavior has been previously observed by G. Mu 
et.al in Ba(Fei_2;Cox)2As2 2j|. A sizeable value of the 
low-temperature specific heat has also been reported 
for Bao.6Ko.4Fe2As2 (70 = 7.7 mJ/mol K^)^ and for 
cuprates superconductors (ssi [36j . 

In general, the origin of the residual 70 observed in su- 
perconducting materials could be caused by pair break- 
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FIG. 1: (Color online) The low-temperature non-lattice part 
of the heat capacity of Ba(Fei_a;Co2:)2As2. The dashed line 
describes the normal state specific heat (see text). Inset: dop- 
ing dependence of the residual specific heat (70), normal state 
specific heat 7n and the dilTerence (771-70). 



ing effects of an unconventional superconductor, crys- 
tallographic defects and disorder, and/or spin glass be- 
havior. It is known that for unconventional supercon- 
ductors the non-magnetic defect and impurities destroy 
the singularity of the gap at nodes due to breaking the 
translational symmetry. This process results in a finite 
density of states induced at the Fermi level. However, 
the presence of the nodal gap imposed by symmetry, like 
in cuprates, is ruled out by the vanishingly small resid- 
ual linear term of the thermal conductivity observed in 
Ba(Fei_a;Coa;)2As2[lll Alternatively, a finite linear 
term in the specific heat is most commonly identified with 
regions of the sample which are non-superconducting. To 
support this simplified notion for these samples, in Fig. 2a 
we plot the non-superconducting fraction of the samples 
determined in three ways. First, if the superconduct- 
ing and non-superconducting regions have similar heat 
capacities then the ratio of 70 /7n will be equal to the 
non-superconducting fraction. Additionally, if the super- 
conducting gap structure is unchanged (which we shall 
demonstrate below) then the condensation energy is sim- 
ply equal to A'^nTc, where ^ is a property of the super- 
conducting gap structure. Using the normal state specific 
heat, we extract the condensation energy for all samples 
by integrating the entropy difference of the normal and 
superconducting state. In the case of Ba(Fei_3;Co3;)2As2 
this approach gives U — 13.15, 1270, 230 and 180 mJ/mol 
respectively for x = 0.045, 0.08, 0.103 and 0.105. If 
a portion of the sample is non-superconducting, then 
U will be reduced from its ideal value. Thus by plot- 
ting \-U / A'^nT^ we obtain another measure of the non- 
superconducting fraction. In this case we choose A = 0.22 
so that the estimate of non-superconducting fraction by 
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this measure is the same at a; = 0.08 as that obtained by 
7o/7n- Finally, susceptibility measurements (see Fig. 2b) 
provide the volume of shielded material, which naively 
should represent the fraction of superconducting mate- 
rial. Consequently, 1 — X(2i<')47r, should again repre- 
sent the non-superconducting fraction. Within the un- 
certainty of demagnetization factors and the unknown 
heat capacity of possible non-superconducting regions, 
all three methods (see Fig. 2a) clearly imply that near 
optimal doping a majority of the sample is supercon- 
ducting ( 85-100%), while proceeding to underdoped or 
overdoped samples a significantly smaller fraction is su- 
perconducting (e.g. < 25% for x — 0.045). Consequently, 
for the remainder of our analysis we will use the hypothe- 
sis that 7o /"fn indeed represents the non-superconducting 
regions of the sample (most probably due to inhomogene- 
ity), and will discuss its origin in more detail at the end 
of the paper. 

Fig. 3 displays the electronic part of the specific heat 
of the superconducting portion of Ba(Fei_^Co^)2As2. It 
is obtained by subtracting the normal state contribution, 
together with a small Schottky contribution below 1 K, 
and normalized by (7„-7o). Interestingly, taking into ac- 
count the superconducting fraction of the specific heat 
only, all the curves collapse below T/Tc = 0.7. To eval- 
uate Tc we have used the entropy balance shown by the 
solid red lined in Fig. 3. As can be seen from the fig- 
ure the width of the transition increases away from opti- 
mal doping. The specific heat jump AC/jnTc = 1.65 for 
X = 0.08 decreases away from optimal doping to 1.5, 1.34 
and 1.05 for x = 0.103, 0.105 and 0.045, respectively [3|. 
This observation could result from a small, but notice- 
able spatial distribution of Tc within the crystals. Pre- 
vious analysis has demonstrated that a single s-wave 
gap cannot reproduce the specific heat data at optimal 
doping[30, 31]. This is illustrated by the single gap fit to 
the optimally doped samplefs^ shown as the dashed line 
in Fig. 3. Clearly there are additional low energy excita- 
tions not captured by a single gap , but can be obtained 
by a two-gap model all dopings|30l IsTj . 

Using a similar approach as presented in RefQ, [13 
we have derived A-f{H) = [C{H) - C(0)]/T for all sam- 
ples studied. It is shown in the inset of Fig. 3. The data 
have been presented in the form A7(iJ)/(7„ — 70) versus 
H/Hc2- The values of Hc2 = 39, 26, 24 and 8.7 T, respec- 
tively for X = 0.08, 0.103, 0.105 and 0.045, have been ob- 
tained from the slope of the upper critical field measured 
by specific heat and the relation Ha2^0.69'^Tc^. It 
should be noted that this relation is only a rough ap- 
proximation strictly valid for a single s-wave gap model. 
However, values obtained are reasonable and agree well 
with those from Refii. As can be seen, all curves show 
roughly the same behavior. It does not match with the 
behavior expected for a simple s-wave order parameter 
where the localized states in vortex cores induce A7 to 
be proportional to H/Hc2- The specific heat is chang- 
ing faster with field than expected for a simple s-wave 
gap scenario. On the other hand, an anisotropic gap will 
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FIG. 2: (Color online) (a) Doping dependence of the non- 
superconducting fraction of Ba(Fei_a;Coa;)2 As2 as determined 
by 7o/7n (circles), 1 — X(2if)47r (squares) and l-U/A'y„T'^ (tri- 
angles), (b) Temperature dependence of magnetic suscepti- 
bility of Ba(Fei_a;Co2,)2As2. Measurements were made in a 
field of 20 Oe applied parallel to the ab plane of the single 
crystals. 



cause the specific heat to deviate from the _ff-linear de- 
pendence. A clean c?-wave superconductor, for example, 
gives A7 cx x/ff/^H- For ah doj )ings we found 
good agreement with a two gap model[4l| with a ratio 
of A„ii„/A,„aa; = 0.5 (see the dashed green line in the 
inset of Fig. 3) which is also consistent with the observed 
temperature dependence. Recently, a similar field depen- 
dence of the electronic specific heat has been obtained by 
Y. Bang within an s± model(42l| with impurity scattering 
and a gap size ratio Asmaii / Aiarge = 0.5 4Jj]. Our spe- 
cific heat does not allow us to resolve whether the smaller 
gap is uniform or possibly even contains accidental nodes. 
Consequently, we cannot comment on the doping evolu- 
tion of small changes to the gap structure which may 
include the lifting of a small nodal component at low 
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energies. However, we can make concrete statements on 
the lack of doping dependence of the major energy scales. 
Within a two gap analysis of our specific heat data the 
ratio of the smaller gap to Tc which controls the low tem- 
perature and low field properties does not vary by more 
than 10% over the doping range studied. The ratio of the 
larger gap to Tc, which controls the size of the specific 
heat jump, may be reduced by as much as 30% in the un- 
derdoped sample and by 15% in the overdoped samples 
both relative to optimal doping. However, we cannot rule 
out the possibility that this small doping dependent evo- 
lution of the larger gap is an artifact created by a spread 
in actual doping concentrations. 

These lack of doping dependence to the gap structure 
is in apparent contradiction with results from Raman 
and thermal conductivity studies on Co doped BaFe2As2, 
which indicate a dramatic evolution of the superconduct- 
ing gap upon doping[2T, ^25^. This may be another mani- 
festation of the extreme sensitivity to sample dependence 
which appears to be common to the Fe-based pnictide 
superconductors, although alternative possibilities exist. 
As alluded to above, the thermal conductivity was mea- 
sured down to 50 mK, a significantly lower energy scale 
than the 0.4 K low temperature limit of our measure- 
ments. Consequently, the apparent discrepancy could 
arise from a change in the gap structure below 0.4 K, 
which we are not sensitive to. Another possibility, which 
was suggested in a recent P-doped BaFe2As2 study that 
encounters a similar contradiction, is that thermal con- 
ductivity is more sensitive to light electron pockets which 
possess the changing near nodal gap structure, while spe- 
cific heat is more sensitive to the fully gapped bands 
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FIG. 3: (Color online) The normalized temperature de- 
pendence of the superconducting state specific heat of 
Ba(Fei_a;Coa;)2As2. The dashed grey line represent the spe- 
cific heat of a single s-wave gap. Inset: field-induced change 
in low temperature specific heat obtained at K for II||c 
by extrapolating the experimental data to zero temperatures 
(see text). The green dashed line is theoretical curve for 
A 

min / ^max — 0.5. 



which possess a higher density of states |43|. Modeling 
within an s± gap structure is also able to reconcile some 
of the apparent discrepancy [l^l. More work is needed to 
determine whether the Raman results can also be under- 
stood in this fashion. Additional doping dependent stud- 
ies by alternative gap sensitive techniques are required to 
help resolve these apparent discrepancies. 

Before concluding, we return to the origin of the resid- 
ual linear term of the specific heat. It is tempting to 
simply attribute the residual linear term in the mea- 
sured crystals as due to "poor" crystals which possess 
phase separation on a length scale smaller than that of 
the EDX probe, and that other measurements were made 
on better crystals. While acknowledging this possibility, 
we note that where comparisons are available between 
these crystals and those measured by other groups grown 
in different laboratories we are in reasonable quantitative 
agreement [23 . Isil |37| . Consequently, we believe the resid- 
ual linear term and its doping dependence are intrinsic 
features of the Co-doped BaFe2As2 system. A residual 
linear term in the specific heat may be a result of gapless 
fermionic excitations or from a distribution of two-level 
systems found in glasses. The magnitude of the residual 
linear term is too large to be accounted for by a structural 
glass [i^l, and a lack of magnetic moments in the over- 
doped samplesli^l rule out a spin glass origin. Thus, the 
natural conclusion is that the system is inhomogeneously 
gapped. The inhomogeneity could exist in real space 
(phase separation) or momentum space (referring to sce- 
narios where portions of the Fermi surface are ungapped) . 
The latter appears to be ruled out by the lack of a lin- 
ear term in the thermal transport in zero field (2l], [S^l as 
well as by the lack of full diamagnetic shielding across the 
doping phase diagram (see Fig. 3). Macroscopic real space 
phase separation on the other hand is ruled out by NMR 
results 15^, _45J , as well as sample uniformity as probed by 
our microprobe analysis. However, the NMR lineshape of 
Ba(Fei_xCoa:)2As2 does broaden increasingly as a func- 
tion of doping j45l|. Similar behavior is observed in hole 
doped cuprates and is believed to be a consequence of 
nanoscaled electronic inhomogeneity Strikingly, in 
cuprates a residual linear term is also observed with a 
qualitatively similar doping dependence and whose ori- 
gin is equally perplexing [36;, i47i], which could indicate 
a common origin. We suggest that the origin of this 
linear term is connected with nanoscaled electronic in- 
homogeneity observed in the cuprates, and suggested to 
be present in the bulk of Ba(Fei_2:Coa;)2As2 by Ning et 
al[45j. STM measurements of Ba(Fei_2:Coa;)2As2 con- 
firm the presence of nanoscale electronic inhomogeneity 
without macroscopic phase separation at the surface 17 1. 
Further work is required to fully understand the origin of 
the residual linear term both in the pnictide and cuprate 
superconductors . 

In summary, using the low-temperature specific heat 
and its magnetic field response, we explore details of 
the superconducting state in Ba(Fei_2;Co2;)2As2 at dif- 
ferent (under, optimal and overdoped) doping regimes. 
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By subtracting the lattice specific heat the temperature 
and magnetic field dependence of the electronic specific 
heat has been studied. The temperature and field depen- 
dencies of the superconducting part of electronic specific 
heat exhibit similar behavior for all doping concentra- 
tions. The temperature variation of C'ei (below Tc) as 
well as its field dependence cannot be described by a 
single isotropic s-wave gap, indicating the presence of 
anisotropic gap structure in the system. Indeed, it has 
been shown recently for optimally Co-doped BaFe2As2 
samples that a minimum of two superconducting gaps are 
necessary to describe the temperature dependence of the 
electronic specific heat [10, IMj . Moreover, the lack of dop- 
ing dependence in Ba(Fei_a;Coa;)2As2 indicates that the 
gap structure does not change significantly as a function 
of doping. The significant residual specific heat observed 



in this system (see also RefH^, HO, HH) is attributed to a 
non-superconducting fraction in the sample and suggests 
that nanoscale inhomogeneity may be an important fac- 
tor in Co-doped BaFe2As2. 



Acknowledgments 

Work at Los Alamos National Laboratory was per- 
formed under the auspices of the U.S. Department of 
Energy, Office of Science and supported in part by the 
Los Alamos LDRD program. Research at Oak Ridge 
National Laboratory is sponsored by the Division of Ma- 
terial Sciences and Engineering Office of Basic Energy 
Sciences. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. 

F. Fang, Nature 453, 761 (2008). 
[3] F-C. Hsu, J-Y. Luo, K-W. Yeh, T-K. Chen, T-W Huang, 

R M. Wu, Y-C, Lee, Y-L. Huang, Y-Y. Chu, D-C. Yan, 

and M-K. Wu, PNAS. 105, 14262 (2008). 
[4] S. Medvedev, T. M. McQueen, L A. Troyan, T. Palasyuk, 

M. L Eremets, R. J. Cava, S. Naghavi, F. Casper, V. 

Ksenofontov, G. Wortmann, C. Felser, Nature Mater. 8, 

630 (2009). 

[5] J. H. Tapp, Z. Tang, B. Lv, K. Sasmal, B. Lorenz, P. 

C. W. Chu, and A. M. Guloy, Phys. Rev. B 78, 060505 
(2009). 

[6] D. R. Parker, M. J. Pitcher, P. J. Baker, L Franke, T. 

Lancaster, S. J. Blundell, and S. J. Clarke, Chem. Com- 

mun. 2189 (2009). 
[7] M. Rotter, M. Tegel, and D. Johrendent, Phys. Rev. Lett. 

101, 107006 (2008). 
[8] A. S. Sefat, R. Jin, M. A. McGuire, B. C. Sales, D. J. 

Singh, and D. Mandrus, Phys. Rev. Lett. 101, 117004 

(2008). 

[9] A. Leithe-Jasper, W. Schnelle, C. Geibel, and H. Rosner, 

Phys. Rev. Lett. 101, 207004 (2008). 
[10] L I Mazin and J. Schmalian, Physica C 469, 614 (2009). 
[11] K. Ishida, Y. Nakai, and H. Hosono, J. Phys. Soc. Jpn. 

78, 062001 (2009). 
[12] K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. 

Kontani, and H. Aoki, Phys. Rev. Lett. 101, 087004 

(2008). 

[13] F. Wang, H. Zhai, Y. Ran, A. Vishwanath, and D-H Lee, 

Phys. Rev. Lett. 102, 047005 (2009). 
[14] S. eraser, T. A. Maier, P. J. Hirschfeld, and D. J. 

Scalapino, New Journal of Physics 11, 025016 (2009). 
[15] M.-H. Julien, H. Mayaffre, M. Horvatic, C. Berthier, X. 

D. Zhang, W. Wu, G. F. Chen, N. L. Wang and J. L. 
Luo, EPL, 87, 37001 (2009). 

[16] K. Terashima, Y. Sekiba, J. H. Bowen, K. Nakayama, T. 
Kawahara, T. Sato, P. Richard, Y.-M. Xu, L. J. Li, G. H. 
Cao, Z.-A. Xu, H. Ding, and T. Takahashi, PNAS. 106, 
7330 (2009). 

[17] Y. Yin, M. Zech, T. L. Williams, X. F. Wang, G. Wu, 



X. H. Chen, and J. E. Hoffman, Phys. Rev. Lett. 101, 
097002 (2009). 

[18] R. T. Gordon, N. Ni, C. Martin, M. A. Tanatar, M. D. 

Vannette, H. Kim, G. D. Samolyuk, J. Schmalian, S. 

Nandi, A. Kreyssig, A. L Goldman, J. Q. Yan, S. L. 

Bud'ko, P. C. Canfield, and R. Prozorov, Phys. Rev. Lett. 

102, 127004 (2009). 
[19] T. J. Williams, A. A. Aczel, E. Baggio-Saitovitch, S. L. 

Bud'ko, P. C. Canfield, J. P. Carlo, T. Goko, J. Munevar, 

N. Ni, Y. J. Uemura, W. Yu, and G. M. Luke, Phys. Rev. 

B 80, 094501 (2009). 
[20] F. Ning, K. Ahilan, T. Imai, A. S. Sefat, R. Jin, M. A. 

McGuire, B. C. Sales, and D. Mandrus, J. Phys. Soc. 

Jpn. 77, 103705 (2008). 
[21] M. A. Tanatar, J. P. Reid, H. Shakeripour, X. G. Luo, 

N. Doiron-Leyraud, N. Ni, S. L. Bud'ko, P. C. Canfield, 

R. Prozorov, L. Taillefer, Phys. Rev. Lett. 104, 067002 

(2010). 

[22] J. K. Dong, S. Y. Zhou, T. Y. Guan, X. Qiu, C. Zhang, 
P. Cheng, L. Fang, H. H. Wen, S. Y. Li, Phys. Rev. B 
81, 094520 (2010). 

[23] Y. Machida, K. Tomokuni, T. Isono, K. Izawa, Y. Naka- 
jima, and T. Tamegai, J. Phys. Soc. Jpn. 78, 073705 

(2009) . 

[24] G. Mu, B. Zeng, P. Cheng, Z. Wang, L. Fang, B. Shen, 
L. Shan, C. Ren, H. Wen, Chin. Phys. Lett. 27, 037402 

(2010) . 

[25] B. Muschler, W. Prestel, R. Hackl, T. P. Devereaux, J. 

G. Analytis, Jiun-Haw Chu, 1. R. Fisher, Phys. Rev. B 

80, 180510(R) (2009). 
[26] L L Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 

Phys. Rev. Lett. 101, 057003 (2008). 
[27] K. Seo, B. A. Bernevig, and J. Hu, Phys. Rev. Lett. 101, 

206404 (2008). 

[28] V. Cvetkovic and Z. Tesanovic, Europhys. Lett. 85, 
37002 (2009). 

[29] K. Hashimoto, M. Yamashita, S. Kasahara, Y. Senshu, N. 

Nakata, S. Tonegawa, K. Ikada, A. Serafin, A. Carring- 

ton, T. Terashima, H. Ikeda, T. Shibauchi, Y. Matsuda, 

larXiv:0907.4399vl 
[30] K. Gofryk, A. S. Sefat, E. D. Bauer, M. A. McGuire, B. 

C. Sales, D. Mandrus, J. D. Thompson, and F. Ronning, 



6 



New J. Phys. 12, 023006 (2010). 
[31] F. Hardy, T. Wolf, R. A. Fisher, R. Eder, P. 

Schweiss, P. Adelmann, H. v. Loehneysen, C. Meingast, 

,arXiv:0910.5006vl 
[32] K. Matan, R. Morinaga, K. lida, and T. J. Sato, Phys. 

Rev. B 79, 054526 (2009). 
[33] J. G. Storey, J. W. Loram, J. R. Cooper, Z. Bukowski, 

J. Karpinski, \arXw-J00lMT4i fl 
[34] G. Mu, H. Luo, Z. Wang, L. Shan, C. Ren, and H-H. 

Wen, Phys. Rev. B 79, 174501 (2009). 
[35] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000). 

[36] N. E. Hussey, Adv. Phys. 51, 1685 (2002). 

[37] S. L. Bud'ko, N. Ni, and P. C. Canfield, Phys. Rev. B 
79, 220516 (2009). 

[38] N. R. Werthamer, E. Helfand, an P. C. Hohenberg, Phys. 
Rev. 147 295 (1966). 

[39] N. Ni, M. E. TiUman, J.-Q. Yan, A. Kracher, S. T. Han- 
nahs, S. L. Bud'ko, and P. Canfield, Phys. Rev. B 78, 



214515 (2008). 
[40] G. E. Volovik, JETP Lett. 58, 469 (1993). 
[41] N. Nakai, P. Miranovic, M. Ichioka, and K. Machida, 

Phys. Re v. B 78 , 214515 (2008). 
[42] Y. Bang. larXi^ 912.5049vl 

[43] J. S. Kim, P. J. Hirschfeld, G. R. Stewart, S. Kasa- 

hara, T. Shibauchi, T. Terashima, and Y. Matsuda, 

|arXELl002.3355vl 
[44] H. V. Lohneysen, Phys. Rep. 79, 161 (1981). 
[45] F. Ning, K. Ahilan, T. Imai, A.S. Sefat, R. Jin, M.A. 

McGuire, B.C. Sales, and D. Mandrus, J. Phys. Soc. Jpn. 

78, 013711 (2009). 
[46] P. M. Singer, T. Imai, F. C. Chou, K. Hirota, M. Takaba, 

T. Kakeshita, H. Eisaki, S. Uchida, Phys. Rev. B 72, 

014537 (2005). 

[47] Y. Wang, J. Yan, L. Shan, H-H Wen, Y. Tanabe, T. 
Adachi, and Y. Koike, Phys. Rev. B 76, 064512 (2007). 



